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High-temperature  steam  electrolysis  (HTSE)  is  a  promising  method  for  highly  efficient  large-scale 
hydrogen  production.  The  HTSE  process  not  only  reduces  the  amount  of  thermodynamic  electrical  energy 
requirement  but  also  decreases  the  polarization  losses,  which  improves  the  overall  efficiency  of  hydrogen 
production. 

In  this  paper,  a  two-dimensional  simulation  method  of  the  efficiency  of  the  HTSE  system  integrated  with 
high-temperature  gas-cooled  nuclear  reactor  (HTGR),  which  changes  two  parameters  simultaneously  in  a 
reasonable  range  while  keeping  one  parameter  constant,  was  presented.  Compared  with  one-dimensional 
analysis  method,  the  effects  of  electrical  efficiency  ( ijei  )>  electrolysis  efficiency  ( r]eSt),  and  thermal  efficiency 
(/7th)  on  overall  efficiency  (/j overall )  were  investigated  more  objectively  and  accurately.  Moreover,  the  critical 
concepts  of  r]es  and  r]th  were  put  forward  originally,  which  were  very  important  to  determine  the  optimum 
electrolysis  voltages  and  operation  temperatures  in  the  actual  HTSE  processes.  The  calculated  critical  value 
of  ijes  was  AG(T)/AH(T)  and  the  actual  rjes  should  be  higher  than  the  theoretically  calculated  one  in  order 
to  maintain  the  high  hydrogen  production  efficiency  of  HTSE  system.  Also,  it  was  very  interesting  to  find 
that  the  critical  rjes  was  the  theoretical  maximum  efficiency  in  SOFC  mode.  Furthermore,  the  critical  value 
of  /7th  was  equal  to  the  value  of  ijeh  which  means  the  overall  efficiency  decreases  with  the  r]es  increasing 
if  the  /7th  in  the  actual  HTSE  process  is  less  than  the  critical  value  of  r]th.  Therefore,  it  is  very  important  to 
control  the  rjth  higher  than  the  critical  value  in  the  actual  HTSE  process  to  get  high  overall  system  efficiency. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

High-temperature  steam  electrolysis  (HTSE)  is  a  promising 
method  for  highly  efficient  large-scale  hydrogen  production. 
Operation  at  high  temperatures  reduces  the  electrical  energy 
requirement  for  the  electrolysis  and  also  increases  the  efficiency 
of  the  power-generating  cycle.  In  addition,  high  temperature  sys¬ 
tems  can  promote  electrode  activity  and  lessen  the  overpotential. 
Therefore,  it  is  possible  to  increase  the  electric  current  density 
and  consequently  decrease  the  polarization  losses  at  high  temper¬ 
atures,  which  improves  the  hydrogen  production  density  and  the 
electrolysis  efficiency.  Thus,  the  HTSE  process  is  advantageous  from 
both  thermodynamic  and  kinetic  standpoints  [1-6]. 

High-temperature  gas-cooled  reactor  (HTGR)  represents 
advanced  nuclear  power  plant  with  efficient,  economic  and  safe 
generation  in  the  future  [7].  The  HTGR  coupled  with  the  HTSE 
system  can  be  expected  to  achieve  very  high  and  satisfactory 
hydrogen  production  efficiency  because  it  can  realize  heat-power 
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coordinated  supply  (the  outlet  temperature  of  HTGR  can  be  up 
to  950  °C).  A  10-MW  high-temperature  gas-cooled  test  reactor, 
the  high-temperature  reactor-10  (HTR-10),  was  operated  in  2000 
at  Tsinghua  University,  Beijing.  This  test  reactor  establishes  an 
experimental  base  for  nuclear  process  heat  application.  The  HTSE 
process  is  currently  being  considered  as  one  of  the  potential  heat 
utilization  systems  for  HTR-10.  But  until  now,  few  studies  have 
been  conducted  on  the  practical  limits  for  the  Coverall  of  the  HTSE 
system. 

In  this  paper,  a  method  of  two-dimensional  simulation  of  the 
efficiency  of  the  HTSE  system,  which  changes  two  parameters 
simultaneously  in  a  reasonable  range,  while  keeping  one  param¬ 
eter  constant,  was  presented.  Through  this  research  method,  the 
effects  of  /7ei,  r]es ,  and  rjth  on  77 overall  were  investigated  comprehen¬ 
sively.  And  the  critical  values  of  rjes  and  77th  were  firstly  put  forward, 
which  were  very  important  to  determine  the  optimum  electrolysis 
voltages  and  operation  temperatures  in  the  actual  HTSE  process. 

2.  Efficiency  modeling 

The  efficiency  model  used  in  this  research  and  the  three  key 
parameters,  electrical  efficiency  (/7ei),  electrolysis  efficiency  (r]es)t 


L.  Mingyi  et  al.  /  Journal  of  Power  Sources  183  (2008)  708-712 


709 


Nomenclature 

E 

Nernst  potential  (V) 

F 

Faraday  constant  (C  mol-1 ) 

AG 

Gibbs  free  energy  change  of  the  reaction  (kj  mol-1 ) 

AH 

enthalpy  change  of  the  reaction  (kj  mol-1 ) 

1 

operating  current  density  (A cm-2) 

Q 

heat  (kj) 

Qoverall 

total  thermal  energy  demand  in  the  electrolysis  pro¬ 

cess 

T 

temperature  (K) 

vth 

thermoneutral  potential  (V) 

Ves 

electrolysis  voltage  (V) 

Greek  symbols 

rje  1 

electrical  efficiency  of  HTSE  system 

rjes 

electrolysis  efficiency  of  HTSE  system 

77th 

thermal  efficiency  of  HTSE  system 

77 overall 

overall  efficiency  of  HTSE  system 

and  thermal  efficiency  (^th)  were  studied  and  derived  intensively 
in  our  previous  work  (see  Ref.  [1  ]). 

The  r]e\  referred  to  the  power-generation  efficiency  of  HTGR, 
which  was  mainly  determined  by  the  outlet  temperature  of  HTGR. 
The  rje\  was  given  by 


(i) 


The  efficiency  model  of  high-temperature  steam  electrolysis 
system  is  shown  as 


Coverall  : 


HHV 


AC{T)/r]elr]es  +  Qxh(T)/r]th  -  (AG(r)/»jes>7th)(l  -  ??e s) 


The  total  differential  equation  of  the  Coverall  is  written  as 

Coverall  =  s  + 


drjel 


where 

Coverall 

drje  1 


dr/e 


HHV  AG(T)/t]e 


dr/th 


(4) 


(5) 


[AG(T)/rle[rles  +  Qth(T)/rlth  -  (AG(T)/>jesi?thXl  -  r/es)]2r/2el 


(6) 


Coverall 

dries 

= _ HHV  AG(T)(l/^e|-l/^th) _ 

[AC(r)/»7el??es  +  Qth(T)/^th  -  (AG(T)/?7es>?th)(l  -  ries)]2rils 

(7) 


Coverall 

dr/th 

= _ HHV[Qth(T)  -  (AG(T)/^es)(l  -  >?es)] _ 

[AG(T)/^elIJes  +  Qth(7)/»)th  -  (AG(T)/lfesl?thXl  -  >)e s)]2^h 

(8) 


where  AG  was  the  generated  electrical  energy  and  Q.  was  the  con¬ 
sumed  thermal  energy  used  for  generating  electricity. 

The  rjes  referred  to  the  overall  efficiency  of  the  HTSE  system, 
including  electrochemical  efficiency  (^eiectrochem).  Faraday  effi¬ 
ciency  ( ?7 Faraday )  and  system  efficiency  (System),  ^eiectrochem  was 
E{T)IVes{i,  T),  where  Ves(i,  T )  was  the  electrolysis  voltage  of  the 
HTSE  system  at  a  given  current  density  i  and  operating  temper¬ 
ature  T,  and  E(T)  was  Nernst  potential  of  the  HTSE  system  at 
temperature  T.  77 system  was  the  parasitic  energy  consumption  of  the 
HTSE  system  resulting  from  parasitic  losses,  such  as  resistance  of 
pipeline,  pumping  work,  AC-DC  conversion,  etc.  The  r]es  was  given 
by 

Pies  —  77eletrochem  ^Faraday  77 system  (2) 

For  the  Faraday  and  77 system  were  generally  very  high,  the  r]es 
was  mostly  determined  by  77electrochem.  And  the  E(T)  remained  con¬ 
stant  at  the  specified  operating  temperature.  Therefore  the  r]es  was 
chiefly  decided  by  electrolysis  voltage  (Ves),  which  was  used  to 
control  the  ries  in  the  real  HTSE  system. 

The  77 th  referred  to  the  thermal  utilization  efficiency  of  the  real 
HTGR-HTSE  system.  It  included  the  thermal  exchange  efficiency 
between  the  HTGR  and  the  HTSE  system,  the  heat  dissipation  of  the 
HTSE  system,  the  heat  consumption  for  preheating  excess  steam, 
hydrogen  and  air  as  well  as  the  waste  heat  recycling.  The  r]th  was 
given  by 

QhTSE  T'inlet  —  ^outlet  /q\ 

77th  =  n f -  tTJ 

kiHTGR  1  inlet 

where  Qhtgr  was  the  thermal  energy  from  the  HTGR  system,  Qhtse 
was  the  consumed  thermal  energy  in  the  HTSE  system,  Tinlet  and 
^outlet  were  the  inlet  and  outlet  temperature  of  the  HTSE  system, 
respectively. 


The  range  of  each  variable  parameter  was  assumed  as  fol¬ 
lows:  (1)  77ei :  40-52%;  (2)  r]es:  50-100%;  (3)  77th:  40-90%;  (4)  T: 
500-1000 °C.  The  two-dimensional  simulation  method  was  real¬ 
ized  in  the  MatLab  software  environment. 

3.  Results  and  discussion 

3  A.  Effect  of  T)ei  OTl  77 overall 

3  A  A.  Effect  of  r\ey  at  various  r)es 

Fig.  1  shows  the  effect  of  77el  on  T7overall  at  various  rjes  while 
keeping  r]th  constant  by  simulation  derivation.  Compared  with  the 
previous  work  [1  ]  (r)es  was  fixed  at  100%),  the  effect  of  77el  on  T7overall 
is  more  complicated  at  various  r]es-  There  is  a  critical  value  of  rjes 


710 


L.  Mingyi  et  al.  /  Journal  of  Power  Sources  183  (2008)  708-712 


Fig.  2.  Effect  of  ije \  on  ^overaii  at  various  r]th  (r]es:  90%,  T :  800  °C). 


Fig.  4.  Effect  of  r]es  on  Coverall  at  various  r]th  (r]e\  \  45%,  T :  800  °C). 


which  will  be  introduced  in  detail  in  the  Section  3.4.1.  When  r]es  is 
in  the  range  from  critical  value  to  100%  and  t/el  (HTGR)  increases 
from  40%  to  52%,  the  effect  of  i/el  on  Coverall  decreases  obviously. 
Also,  the  downward  trend  of  r]e\  effect  increases.  On  the  contrary, 
when  the  value  of  ijes  is  in  the  range  from  50%  to  critical  value  and 
i)ei  (HTGR)  increases  from  40%  to  52%,  the  effect  of  rje\  on  Coverall 
increases  but  the  upward  trend  of  i/el  effect  decreases. 

3.1.2.  Effect  of  r]ei  at  various  ijth 

Figs.  2  and  3  show  the  effect  of  rjei  on  Coverall  at  various  rjth  while 
keeping  r]es  constant.  The  effect  of  r]e\  on  Coverall  decreases  slightly 
with  the  increase  of  rje\  from  40%  to  52%  at  various  r]th  from  40%  to 
90%,  when  rje s  is  fixed  above  the  critical  value  (e.g.  90%),  as  shown 
in  Fig.  2.  However,  the  effect  of  r]e\  on  Coverall  increases  with  the 
increase  of  77el  at  various  rjth ,  when  rjes  is  fixed  below  the  critical 
value  (e.g.  70%),  as  shown  in  Fig.  3.  Although  ijei  has  an  opposite 
effect  on  Coverall  at  above  and  below  the  critical  value  of  i)e s,  the 
values  of  r]e\  effect  both  change  slightly. 

3.2.  Effect  of  ties  on  r]ovemii 
3.2.1.  Effect  of  r]es  at  various  r\th 

Fig.  4  shows  the  effect  of  r\^s  on  Coverall  at  various  r]th  while 
keeping  ^el  constant.  The  effect  of  r]es  is  less  significant  than  that 
of  rje\.  When  the  value  of  r]th  is  less  than  that  of  rje\,  the  effect  of 
Ties  on  Coverall  is  negative,  which  indicates  Coverall  decreases  with 
the  increase  of  r]es,  as  shown  in  Fig.  4.  The  same  result  can  also 


be  obtained  from  Eq.  (7).  The  value  of  rjth ,  which  is  equal  to  rjeh 
is  called  the  critical  value  of  r]th.  Therefore,  ijth  should  be  ensured 
higher  than  the  critical  value  of  r]th  in  the  actual  HTSE  process  to 
guarantee  the  high  efficiency  of  hydrogen  production.  In  addition, 
the  effects  of  rjes  on  Coverall  decrease  slightly  with  the  increase  of 
ries  from  50%  to  100%  at  various  r]th  from  the  critical  value  to  90%. 

3.2.2.  Effect  ofries  at  various  rjei 

Fig.  5  shows  the  effect  of  rjes  on  Coverall  at  various  rjei  while  keep¬ 
ing  r]th  constant.  The  effect  of  r]es  on  Coverall  decreases  obviously 
with  the  increase  of  r]es  from  50%  to  100%  at  various  ^el  from  40%  to 
52%.  The  change  trends  of  the  effects  of  rjes  are  similar  at  various  rje\. 
The  change  values  of  rje s  effect  increase  slightly  with  the  increase 
of  T]e 

3.3.  Effect  of  rjth  on  r]overall 
3.3.1.  Effect  ofrjth  at  various  r]es 

Fig.  6  shows  the  effect  of  r]th  on  Coverall  at  various  r]es  while  keep¬ 
ing  r)e\  constant.  r/th  has  the  smallest  effect  among  the  three  key 
parameters  apparently  which  is  consistent  with  the  result  of  the 
previous  work.  As  shown  in  Fig.  6,  there  is  another  critical  value  of 
ries.  When  the  value  of  r]es  is  less  than  the  critical  value,  the  effect 
of  77th  on  77overan  is  negative,  which  means  77overaii  decreases  with 
the  increase  of  77th.  The  critical  value  of  r]es  can  be  determined  to 
be  AG(T)/AH(T)  by  derivation  of  Eq.  (8).  The  derivation  process  is 
described  in  the  Appendix  A.  Therefore,  rjes  should  also  be  ensured 
higher  than  the  critical  value  in  the  actual  HTSE  process  to  guar- 
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Fig.  3.  Effect  of  T]q\  On  Coverall  Various  T^th  (^7es*  70%,  T.  800  C). 


Fig.  5.  Effect  of  rjes  on  Coverall  at  various  rje i  90%,  T :  800°C). 
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antee  the  high  system  efficiency.  In  addition,  the  effects  of  rjth  on 
Coverall  decrease  slightly  at  various  r]e s  from  the  critical  value  to 
100%. 

3.3.2.  Effect  ofr]th  at  various  r]d 

Fig.  7  shows  the  effect  of  rjth  on  Coverall  at  various  ije\  while  keep¬ 
ing  rje s  constant.  The  effect  of  rjth  on  Coverall  decreases  obviously 
with  the  increase  of  rjth  from  40%  to  90%  at  various  rje\  from  40%  to 
52%.  The  change  trends  of  the  effects  of  rjth  are  similar  at  various 
rje\.  The  change  values  of  the  effects  of  rjth  increase  slightly  with  the 
increase  of  ije\.  The  results  of  the  effects  of  rjth  are  similar  to  that  of 
the  effects  of  r)es  at  various  r]e\. 

3.4.  Critical  values  analysis  ofHTSE  system 

At  a  specified  operating  temperature,  rjeh  r)es,  and  rjth  had  signif¬ 
icant  influence  on  Coverall  °f  the  real  HTGR-HTSE  system.  i/el  was  the 
power-generation  efficiency  of  HTGR,  which  was  kept  almost  con¬ 
stant  to  the  specified  HTGR  system.  The  i/el  of  HTR-10  of  Tsinghua 
University  is  about  48%  at  850  °C  (the  outlet  temperature  of  HTGR). 
r)es  was  the  efficiency  of  the  HTSE  system,  which  was  mainly  deter¬ 
mined  by  electrolysis  voltage.  The  electrolysis  voltage  can  be  used 
to  adjust  the  r)es  in  the  real  HTSE  system.  Therefore,  the  electroly¬ 
sis  voltage  is  the  key  parameter  in  the  real  HTSE  system.  rjth  was 
the  thermal  utilization  efficiency  of  the  real  HTGR-HTSE  system 
determined  by  the  outlet  temperatures  of  HTGR  and  HTSE  systems. 
77th  increased  with  the  increase  of  the  outlet  temperature  of  HTGR 
system  and  the  decrease  of  the  outlet  temperature  ofHTSE  system. 


Fig.  7.  Effect  of  rj^Yi  Oil  Coverall  VcirioUS  T]q\  ( rj eg.  90/^,  T.  800  C). 


Fig.  8.  Critical  values  of  r]es  at  various  temperatures. 

Generally,  77th  was  more  than  80%  in  the  real  systems  because  the 
waste  heat  could  be  recycled,  such  as  vaporizing  water,  preheating 
the  steam  and  hydrogen  and  so  on. 

There  were  two  critical  values  of  r]es  and  77th  in  the  real 
HTGR-HTSE  system,  as  was  mentioned  above.  These  two  critical 
values  were  very  important  to  determine  the  optimum  process 
parameters  in  the  actual  HTGR-HTSE  process.  The  critical  value  of 
ties  was  the  minimum  electrolysis  efficiency  in  the  real  HTGR-HTSE 
process,  which  was  controlled  by  Ves-  Through  the  critical  value 
of  rjes,  the  critical  value  of  Ves,  which  is  a  key  parameter  in  the 
actual  HTSE  operation,  could  be  calculated.  The  critical  value  of  Ves* 
which  could  be  further  recommended  as  the  maximum  electrolysis 
voltage  in  the  actual  HTGR-HTSE  process  to  ensure  highly  efficient 
operation,  was  an  equilibrium  point  of  the  hydrogen  production 
density  and  electrolysis  efficiency  of  the  system.  With  the  increase 
of  Ves,  the  hydrogen  production  density  increases  while  the  elec¬ 
trolysis  efficiency  decreases.  Conversely,  it  is  the  case  similarly.  The 
critical  value  of  77th  was  the  minimum  thermal  efficiency  in  the  real 
HTGR-HTSE  process,  which  was  mainly  adjusted  through  the  outlet 
temperature  of  HTSE  system. 

3.4.1.  Critical  values  ofrjes  and  Ves 

At  the  critical  value  of  r)es,  the  effects  of  77el  on  77overall 
( Coverall /d^ei)  kept  constant  with  rje\  changing,  which  meant  the 
Coverall  remained  unchanged  with  77el  increasing,  as  shown  in  Sec¬ 
tion  3.1.1.  Then  the  critical  value  of  ijes  could  be  derived  by  Eq. 
(6).  For  Coverall /^4ei  kept  constant,  the  second-order  differential 
equation  (92 Coverall /d^)  was  zero,  the  critical  value  of  r]es,  which 
was  AG(T)/AH(T),  could  be  deduced.  The  derivation  process  was 
described  in  the  Appendix  B. 

At  the  critical  value  of  r]eSt  the  effect  of  ij th  on  77overan 
( Coverall /d^th)  was  zero,  which  meant  the  Coverall  decreased  with 
r]th  increasing,  when  the  f]es  in  the  actual  HTSE  process  was  less 
than  the  critical  value  of  r]es,  as  shown  in  Section  3.3.1.  Therefore, 
the  r)es  in  the  actual  HTSE  process  should  be  controlled  to  be  higher 
than  the  critical  value  to  maintain  the  high  efficiency  of  system. 

As  A G(T)  and  A H(T)  were  both  the  functions  of  temperatures, 
the  critical  value  of  r]es  should  also  be  the  functions  of  tempera¬ 
tures.  Fig.  8  showed  the  critical  value  of  rjes  at  various  temperatures. 
There  was  an  interesting  result  that  the  critical  value  of  r]es  was  the 
theoretical  maximum  efficiency  in  fuel  cell  mode. 

r)es  was  controlled  by  Ves,  which  was  the  actual  operational 
parameter  in  the  real  HTSE  process.  The  critical  value  of  Ves  could 
be  calculated  by  the  critical  value  of  r]es  through  the  equation  of 
rjes=E(T)lVes(T).  Fig.  9  showed  the  critical  values  of  Ves  at  various 
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temperatures.  The  critical  values  of  Ves  were  the  recommended 
maximum  electrolysis  voltage  in  the  real  HTSE  process,  which 
increased  from  1.275  to  1.293  V  with  temperatures  increasing  from 
500  to  1000  °C. 

3.4.2.  Critical  value  of  rjth 

At  the  critical  value  of  rjth,  the  effects  of  rjes  on  Coverall 
(d?7 overall /d?7es)  was  zero,  which  meant  the  Coverall  decreased  with 
the  rje s  increasing,  when  the  rjth  in  the  actual  HTSE  process  was  less 
than  the  critical  value  of  rjth ,  as  shown  in  Section  3.2.1.  Therefore, 
the  rjth  in  the  actual  HTSE  process  should  be  controlled  to  be  higher 
than  the  critical  value  to  maintain  the  high  efficiency  of  hydrogen 
production.  Otherwise,  the  overall  efficiency  of  HTGR-HTSE  sys¬ 
tem  would  decrease  with  the  r]es  increasing.  The  critical  value  of 
r]th  could  be  derived  from  Eq.  (7),  which  was  equal  to  the  value  of 
r]e\.  Comparison  to  rjeh  which  was  less  than  50%,  rjth ,  which  was 
more  than  80%  in  the  real  systems,  was  much  higher.  Therefore,  the 
rjth  was  higher  than  the  critical  value  generally  in  the  actual  HTSE 
process. 

4.  Conclusions 

In  this  paper,  a  two-dimensional  simulation  method  of  the  effi¬ 
ciency  of  the  HTGR  coupled  with  HTSE  system  was  presented.  The 
effects  of  electrical  efficiency  (?7el),  electrolysis  efficiency  [r]e s,),  and 
thermal  efficiency  (r/th)  on  overall  efficiency  (^overall) were  investi¬ 
gated  more  objectively  and  accurately.  And  the  determining  factors 
of  these  three  different  efficiencies  in  the  real  HTGR-HTSE  system 
were  also  introduced. 

The  critical  concepts  of  r]es  and  r)th  were  put  forward  origi¬ 
nally  in  the  real  HTGR-HTSE  system,  which  were  very  important  to 
determine  the  optimum  electrolysis  voltages  and  operation  tem¬ 
peratures  during  the  actual  HTSE  processes. 

The  calculated  critical  value  of  r]es  was  AG(T)/AH(T)  and  the 
actual  r]es  should  be  higher  than  the  theoretically  calculated  one 
in  order  to  guarantee  the  high  hydrogen  production  efficiency  of 
HTSE  system.  Also,  it  was  very  interesting  to  find  that  the  critical  r]es 
was  the  theoretical  maximum  efficiency  in  SOFC  mode.  The  critical 
value  of  Ves*  which  was  determined  by  the  critical  value  of  r]es,  was 
a  key  operational  parameter  in  the  real  HTGR-HTSE  system.  It  was 
recommended  as  the  maximum  electrolysis  voltage  to  ensure  high 
hydrogen  production  efficiency  during  the  actual  HTSE  processes. 


The  critical  value  of  rjth  was  equal  to  the  value  of  r]e\.  Generally, 
the  actual  r]th  was  much  higher  than  rje\.  Therefore,  the  rjth  was 
higher  than  the  critical  value  in  the  real  HTGR-HTSE  system. 
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Appendix  A.  The  derivation  process  of  the  critical  value  of 
rjes  related  to  rjth 

For  Coverall /d^th  =  0  at  the  critical  value  of  r]e s,  then  from  Eq.  (8) 
Qth(T)  -  (AC(T)/»7es)(l  -  rjes)  =  0  (Al) 

Then 

AG(T)/ r)es  =  A  G(T)  +  Qth(T)  =  A  H(T)  (A2) 

So  the  critical  value  of  r]es  =  AG(T)/AH(T). 

Appendix  B.  The  derivation  process  of  the  critical  value  of 
rjes  related  to  rjei 

For  the  effect  of  rje\  on  Coverall  is  constant  at  this  critical  value, 
S^overall/^e]  is  zera 

From  Eq.  (6),  the  second-order  differential  equation  of  the  Coverall 
is  given  by 

9* 2 3 4 5 6 7>? overall  -2  HHVAG(T)  A  AG(p\ 

9>?e  l  *?es7?g|t2  v  V  eldest) 

where 

t=Acm  +  Q*(I)_AC(r)(1  es) 

^el^es  r]th  T]esr]th 

For  -2  HHV  AG{T)/r]esr]^t2  ±  0,  then  we  can  obtain 

,  Acm_0 

*1eMest 

Substituting  t  into  Eq.  (B2),  we  can  obtain 
A  G(T) 

ACCT)  +  QthCn»?elW>?th  -  AC(T)/?el(l  -  VesVrith 
Then,  we  can  derive 

QthCn»?elW>?th  -  AG(T)'?el(l  -  Ves)/r]th  =  0  (B4) 

So  the  critical  value  of  r]es  =  AG(r)/AH(T). 
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